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SUBJECT: BIWEEKLY REPORT, OCTOBER 19 ^ 1953 

To: Jay Wo Forrester 

From: Scientific and Engineering -Computation Group 

1<. MATHEMATICS. CODING AND APPLICATIONS 
lol Introduction 



During the period covered by this report 207 coded programs were run 
on the time allocated to the Scientific and Engineering Computation (S&EC) Group, 
These programs represent part of the work that has been carried on in 20 of the 
problems that have been accepted by the S&EC Group o Progress on each of these 
problems is given below in terms of programming hours g minutes of computer time , 
and progress reports as submitted by the programmers in question o 

The following five problems have been accepted by the S&EC Group for 
solution on WWIi 



Problem Noo lUB 

Problem Noo 149 

Problem Noo 151 

Problem Noo 152 

Problem Noo 1$3 



Elliptic Boundary ¥alue Problem (Lifting Line Theory) 
Ho Glantzg, Mathematics Department 

Digital Methods of Detecting Signal From Noise 
G . Dinneen ^ Lincoln Laboratory 

Programming and Analysis of General Game of NIH 
Do Sternlightg Mathematics Department 

Diffusion In An Oxide Coated Cathode 

Bo Frost, Digital Computer Laboratory 

Gust Responses Simultaneous Linear Integro'iiMff <i*^lquatiibn 
Ko Fossg Aeroelastic and Structures Laboratory 



Detailed reports describing these problems will be presented in future biweekly 
reports o 

Although none of the programmers associated with the S&EC Group is developing 
routines expressly for the library of subroutines j some of the programs that ha^ve 
been developed for specific applications have proved to be of general utility <, 
Such programs are transcribed to a subroutine form^ tested, and then written up for 
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the library o At present 13 such routines are available with two more \mder test* 
Progress in this work is described under Problem Noo 141« This library does not 
include the automatic (output and PA) routines in the CS« 



1«2 Programs and Computer Operation 
100, Go 



Comprehensive System of Service Routines; Gombelic, 1»5 hours} Benman, 
59 hour^; Demurjian, S. 5 hours? Frankovich, 8o5 hours; Kopley, 6 hours; 
Porter, 4 hours; Helwig^^lV'hcturS; WWI, 224 minutes 

Tests were made for using the character generator in CS scope post mortemso 
Results revealed that this method is about six and a E&lf times as fast as the 
previously used point-by-point method*, All 2016 registers of storage were recorded 
as octal constants in 11*5 seconds o However, it has not yet been decided whether 
the numbers are sufficiently legible o The size of the numbers was about twice that 
used by Group 61 o This was accomplished by adjusting thegaino Whether this cm 
be done without affecting point-by-point plotting is being investigated » 

Kopley 

The new conversionrroutine for transforming (24.>6) floating-point numbers 
into normalized decimal form (+ ol234-5678 ( + 12), written for the new drim Post 
Mortem, was completed, and the first tests indicate that it operated correctly » 
It is written in closed subroutine form, in WWI code, and is entered with mag- 
netic-tape imit #3 running and with ca x in the AC, where x is the address of 
the first half of the number in storage to be conyertedo The routine now occupies 
about 200 registers (including its own programmed arithmetic) and requires about 
116 milliseconds to record the 15 Flexowriter characters per niimber on tape for 
delayed printing. No round-off facility has been included o 

This routine is to be rewritten as a sub-routine for users of GS to print 
out the contents of the MElAo 

Denman 

It was discovered that some of the numbers printed out by the automatic GS 
output routines were in error in the fifth decimal place* ^his error was due to 
a special-add overflow that was not picked upo The error has been correctedo 

. Porter 

101 o Optical Properties of Thin Metal Films ; Denman, 2 hours; Loeb, 3 hours; 
Richmond, 4 hours; WWI, 7 minutes 

Although the program worked with one set of parameters, it failed to work 
with another set of parameters of the same order of magnitude » There was a divide 
and an 'its alarm. We are now trying to find the source of the trouble e 

Loeb 

106o MIT Seismic Project; Briscoe, 4.0 hours; Simpson, 20 hours; Walsh, 20 hours; 
WWI, 142 minutes 

Correlations have been computed for intervals of seismogram over which linear 
operators have been computed o These can be interpreted so as to optimize the 
operator structure before least squares fitting is effectedo A program to set up 
normal equations for general operator structure has been written and is being 
tested o 
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Another study underway is the comparison of averaged operators and averaged 
error curves o New programs include a subroutine to plot numbers as well as print 
them out on the delayed printer, and an averaging program to sample the changing 
mean on seismogram readings.. 

Further tests have been made on a power-density-spectrum programo 

Robinson 

108 o An Interpretive Program ; Hazel, 1 hour; Zierler, 60 hours; WWI, 35 minutes 

Enough of the program now operates satisfactorily so that it can be used for 
the solution of non-trivial problems o The following are the only features not yet 
operating correctly: automatic computation of special functions and ihe assignment 
of variables to the drum. It is expected that the entire program will be in oper- 
ation and that a report describing the program in detail will appear within the 
next few weeks. 

Zierler 

109 o li ghter Gunsieht Calibration. 8th Order DoEo ! Hellman, 1 hour; Frankovich, 
o5 hours; WWI, 13 minutes 

Three runs were made during this periodo An error was encoimtered in the 
CS output routine. This made it necessary to change the program print routine to 
incorporate Tape No. 2299-6. A test run incorporating this print subroutine was 
successful. 

^he initial conditions for a three-dimensional pursuit-course problem are 
being determined. Coding for this problem should be completed during the next 
biweekly period. 

Hellman 

116. Torpedo Impulse Response; Convolution ; Hamilton, 30 hours; WWI, 11 minutes 

All convolutions and computations have been completed for the report being 
prepared on the problem of finding impulse responses for various input-output pairs. 

An inverse transform and possibly another convolution will be performed on 
WWI after the report has been completed. After these have been performed ouf 
problem will be terminated. 

Kramer 

119 o Spherical Wave Propagation; Ralston, 20 hours; WWI^ 102 minutes 

One successful run of 37 minutes has carried the computation twice ^as far 
as previously. A first glance at the results tends to indicate that a shock is 
not going to form as expected. 

The next run will continue the computation from where it stopped and may 
indicate conclusively whether or not a shock will form. If so, this will complete 
this phase of the problem. 

Ralston 

126. Data Reduction; ^oss, 60 hours; Hamilton, 10 hours | Frankovich, 2 hour®| 
WWI, 193 minutes 
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Several sections of the Mistake Diagnosis Routine are now known to be 
operative o Considerable time was lost during this period due to the difficulty 
of locating a tape-preparation mistake » Further time was lost due to the computer 
dropping a sign digit repeatedlyo For this no explanation has been fotrnd* Several 
necessary modifications to the program disclosed by a successful test run have not 
been tested because of tape errors. 

Ross 

132, Subroutine Study for the Numerically Controlled Milling Machine; Runyon, 
35 hours; Frankovich, 4. hours? WWI, lO^L minutes 

Writing of the program for preparing NCMM tape for a super-sonic nozzle was 
completed* Another program which furnishes tool-center coordinates of one cross 
section of the nozzle was rim successfully using the nozzle dimensions o Two sub- 
routines, one of which is used in the iiozzle program, are in the process of being 
tested© 

R^anyon 

1360 Matrix Equations ; Arden, 20 hoiirs; Helvi^il6 hours | WWI, 21 minutes 

The program described in the biweekly report of June 1, 1953 has been 
transcribed into subroutine form. After ftirther testing, the transcription will 
be placed in the subroutine library of WWIo 

Arden 

An n step procedure for solving systems of linear equations has been derived 
which includes the conjugate gradient method and the method described by Craig 
(M-2229) as special cases* The generalization suggests other n step methods but 
reasons for using , these methods in place of the conjugate gradient method have 
not been found© work on further generalization is continuing o 

Helwig 

137 o Investigation of Atmospheric Turbulences Autocorrelation. Crosscorr elation 
and Fotirier Transforms ; Summers, 20 hours | Block 8 hours; Kopley, o5 hours; 
WWI, 116 minutes 

Successful determination of the Run II Autocorrelation fimctions for lOO-- 
200 values of the original data time spacing (4^ t) was accomplished using the 
basic crosscorrelation program (T-2751 or 3068), but starting the second data- 
tape read-in on the fifth block© For R\m I, however, the 100-200 At points were 
not consistent with the 0-100 At computation^ This may have been due to startiig 
with the wrong block; the proper starting points will be clearly marked on the 
tapes and performance requests resubmitted© 

The Fourier transform program CT-2235ib18) was used to take the transform 
of a unit step and results seemed satisfactory© fhis program was then modified 
to eliminate the use of Simpson's rule, but, when used on the flight data, *»neg- 
ative" power was still obtained© The program will have to be examined more 
closely© 

ompletion of the problem now depends almost entirely on suitable modification 
of the Fourier transform program so as to eliminate that *^negative" power© 

Summers 
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138 <, Spheroidal Wave Fimctions ; OiombeliCg 4^5 hours? Little, 30 hours | Corbatij, 
50 hours I WI, 84 minutes 

The "Ratio Obnverters and Accumulators" program has been testedo A run of 
this routine with the previously tested "Layout" program utilizing the magnetic 
drum has now been made and works satisfactorily o 

A (15,0) decimal integer delayed print routine was required for testing. 
3uch a routine was assembled, tested and submitted to the library of subroutines 
for general use* ^ 

Gorbato 

The routine for calculation of the ratios of the coefficients is being 
modified to give greater accuracy and to take less time* It is about ready for 
useo In the course of testing an error was foimd in the auxiliary buffer subroutine, 
File 2767; this mistake has been corrected <» 

Little 

140, Summer Session System: Combelic, 48 hours | Helwig, 8 hours 5 Franko^ich, 
20 o 5 hours; Hoy, 8 hours; Siegel, 46 hours; WWI, 335 minutes 

The changes in the input-output instructions of the Summer Session Computer, 
which were described in the previous biweekly report, have been completed and 
testedo These changes are now in effects 

A test program, which tests every instruction of the computer ^nd indicates 
which, if any, of them is not operating properly has been written* Some minor 
errors, located by this program and hy other means, have been corrected* Effort 
is being directed toward locating and eliminating all such remaining errors » 

The several punched tapes which comprise the computer have been combined 
into one tapeo fhis permits reading the tape into HWI without special instruct- 
ions to the operators* A program which will record the subroutines on magnetic 
tape imit 1 and will permit their automatic reading from the magnetic tape to the 
drum is being written* , ^his program will obviate entirely the need for reading 
in a long punched tape each time the computer is to be usedo 

Thought is being given to the development of a subroutine library for the 
Summer Session Computer* Since there is no provision for the use of relative 
addresses, means for incorporating library subroutines into a program without 
conflict of floating addresses are being studied* It has been decided to reserve 
floating addresses numbered 900 or greater for the use of library subroutines* 

A summer session subroutine for finding the square root of a number is being 
written and will undergo tests dtaring the ne:^t biweekly period* Another subroutine 
written by C* W.Adams, for alphanumerical character display on the oscilloscope has 
been tested and operates properly* This routine will be available for use by 
students in Course 6*535* 

A pfogram, also the work of Professor Adams, for the solution by several 
methods of the differential equation y' « 1 - y^ was run on the computer with 
good results* 

A study of the interpretive subroutines used by the computer with the aim 
of producing accurate annotated copies of their contents has been initiated* No 
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such copies exist at presento They will be an invaluable aid in amending or 
modifying the operation of the computer in the future if changes shoxild appear 
advisable o 

Siege 1 

141 o S&EC Subroutine Study ; Arden, 4- hours j Vonderburgh, 2 hours 5 WWI, 7 minutes 

A new subroutine, consisting of a 23-register addition to subroutine tape 
2756m5 is being testedo This addition will enable the user to arrange the output 
of 2756m5 in certain simple formats* 

Arden 

A new routine has been added to the subroutine library bringing the total 
to 13 o The new routine will calculate the eigenvalue with the largest absolute 
magnitude of a real matrix all of whose eigenvalues must be realo This routine 
was developed under Problem No<, 146» 

A subroutine for solving a system of linear algebraic equations (the number 
of equations is limited physically by the available magnetic drum space and 
mathematically by the accumulated roundoff errors) is being tested by Do Arden 
(Problem No* 136)0 A subroutine for the delayed (ioeoj via magnetic tape) print 
out of (15 fO) decimal integers is being tested by Fo Corbato (Problem No« 138) o 

Vanderburgh 

1/42 o A Studt of Shock Waves ; Sydney j 60 hours | Bart, 60 hours | WWI, 4-77 minutes 

^n the past two weeks the program for the two-dimensional case has been 
unde r continued te sting « 

The program was run in its entirety for the first time on October 18 5 195 3> 
and appeared to be functioning satisfactorily o The results of the run are being 
checked and if they are correct, the program will be ready for further production 
runs* 

Sydney 

14.30 Vibrational Frequency Spectrum of a Copper Crystal g Third Order Polvnomial g 
Jacobsen, 12 hours | WI9 17 minutes 

The complete program in interpretive (24.96) is nearly finished and will be 
tested with a coarse mesh of 95 points o However, I believe a finer mesh of 
25,4-95 points will be required ultimately o %is being the case 9 I shall have to 
re-program the routine in PA (15,15) in order to save coB^uter machine timeo 
The object of the interpretive (24.96) routine above was to establish a standard 
of accuracy and to learn the nature of 'lay problem in terms of Whirlwind language « 

Jacobsen 

l^Vo Self -consistent Molecular Orbitals ; Arden, 7o5 hours g^Meckler, 15 hours g 
WWI?s> 56 minutes 

The program described in the biweekly report of September 20 j 1953 has been 
written o Several technical errors have been discovered and corrected e 

Heckler 
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1460 Largest Eigenvalxae of Realg Symmetric Hatrlx s Siegelg 1 hoiir^ Temkinj 2 
hours I Wis, 14 minutes 

The problem has now been completed and the subroutine for finding the 
largest eigenvalue of a real matrix (all of whose eigenvalues are real) entered 
into the library of subroutines as #3l68ml2e The program was checked against a 
known problem taken from Milne (Numerical Calculus) and was correct to as many 
digits as taken there., Another problem was checked against Meckler's program 
(Problem #134» tape 26l3p2m3) and the results checked to six significant figures, 
and differed by one in the seventh place « 

Temkin 

147 o Energy Bands in Crystals ; Arden, I4 hours § Howarth? 20 hours | WWIj, 55 minutes 

A program for the solution of the differential equation and the calculation 
of the expressions v^^ (E,r) described in the biweekly report of October 5> 1953> 
has been successfully completedo A program for the calculation of the F(E9Eo»r) 
has also been written and is being testedo 

Arden 

lo3 Operating Statistics 
Computer Time 

The following indicates the distribution of WWI time allocated to the S&E6 

Group o 

Programs 33 hours 9 48 minutes 

Conversion 11 hours ^ 44 minutes 

Magnetic-dr\3m Test 05 minutes 

Scope Calibration ,__ 62 minutes 

^otal Time Used 46 hours ^ 39 minutes 

Total Time Assigned 55 hours ^ 52 minutes 

Usable time, percentage 83e56^ 

Nixmber of Programs 207' 

1<»4 Summary of Tape Room Bulletin Board Memoranda (lo Hazel) 

(These memos are intended to inform programmers of changes in coding 
procedure, WWI facilities, etco) 

Programmed Arithmetic (GS) 

An error which is not easily correctable has been discovered in the Programmed 
Arithmetic for the instruction ics b e Program^lers using the buffer or the Auxiliary 
Buffer Subroutine, #2767, should refrain^ until further notice, from using the ics 
instruction on these buffer registers <, Suggested alternatives for ics b ares 
lo ica (RC + „0) 2e ica b 

ii^ b ^ ^^^ ""^'^^ 

i dv Su b-Elbck of the P_A. 

The i dv sub--block''of the PA will be punched out on the 556 tape every time 
the PA is used whether or not an i dv occurs in the program* 
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2o COMPUTER ENGINEERING 

2«1 WWI System Operation 

2 oil Core Memory (No Lo Daggett) 

On October 6, the Core Memory ended a six-day stretch without any parity 
alarms Since that time there have been several failures in Bank Bo Bank Aj, 
the new core bank, has still given no failures since the 6th o It appears that 
the higher delta noise encountered with the old (Bank B) cores makes this bank 
more susceptible to power transients© 

(Lo Lo Holmes 9 AoJ«Roberts) 

We have continued to remove or reroute the power wiring and video cabling 
in the EXl ~ EX8 areao The work was completed on October 10 o The racks and 
wire -ways in this region will be dismantled and removed dxiring the forthcoming 
biweekly period o This will allow us to locate Core Memory in that section of 
the roomo 

2«12 Auxilikry Magnetic--Drum System (KoEoMcVicar) 

The problem of slupeish relay operation in the group-selection system has been 
attacked via two routes© The first (and so far not very successful) alternative 
was to attempt to increase the relay operate speed o The second approach has been 
to increase the relay setup time counted by the In-Out Delay Counter o 

There is reason to suspect that the auxiliary drim is occasionally making 
an error in a single digit and a single register© The error is permanent rather 
than transient, ioS", the drum surface has been altered o The digit and register 
In error seem to vary and are not confined to any particular group o 

Tests on the auxiliary drum with an interlace of 4 ha¥e been run^ and marginal 
checking of the drum circuits indicates that satisfactory operation can be expect- 
ed © Permanent change of , the interlace will await elimination of the errors ment- 
ioned above to avoid confusion which might result from piliiag the troubles one 
upon the other o 

Work on the logic for a drum parity check is in progress o In the meantime 
we are starting construction of the elements in the driam required for a parity 
system© 

2 ©2 Terminal Equipment 

2 ©21 Ferranti PETE (FoE© Irish) 

A clutch-brake driver for the Ferranti PETR has been constructed© Using this 
circuit along with a test-equipment setup 3 the new PETR has been operated in what 
is essentially a line-by-line mode© 
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The maximum speed of the paper t ape when used in this reader is 25 inches of 
tape or 250 Flexowriter characters per second o If the brake is applied while 
the unit is operating at this speed, the tape will stop its motion in from Oo5 
to loO milliseconds o Starting takes from Oo75 to lo25 millisecondso 

A design for the information-channel amplifiers to be used in this reader 
is being developedo Houghly speaking, each channel aniplifier will consist of 
a phototube, followed by a single triode amplifier stage and a cathode-follower 
output o 

2o22 Magnetic Tape (EoPoFarnsworth) 

Minnesota Mining Coo has supplied us with several reels of high-output green 
tape on the Scotch acetate base for trial o We have installed one of these reels 
on unit 3B to check for the twisting action which developed in the first product- 
ion batch. This tape was brcyken, apparently through mishandling, less than 12 
hours after installation, reempha sizing the importance of the nQrlar baseo 3M is 
still trying to get lo5-inil n^rlar in 1200-foot lengths from Duponto 

Inuring this period, three Raytheon read/record heads become marginal due to 
mechanical -wear, and one failed due to insulation breakdown o Replacement of these 
heads necessitated rerecording of the "permanently" recorded data on units and 
1, much to the consternation of individuals who had assumed greater permanence in 
magnetic tape than is attributed to punched paper tape? in spite of the possibili» 
ties of tape damage due to mishandling, accidental erasure (except on unit which 
is locked in "read"), head failures, etco 
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4.0 ADMINISTRATION AND PERSOMML 

Terminated Staff (Jo Co Proctor) 
Roy Eo Hegler 

New Non-Staff (R« A« Osborne) 

Elinor Albanese is a new messenger at the Barta Building o 

Pauline Gochary is Mro Morley's new secretary o 

Robert Johnston has joined the ■Whittemore Building janitor crewo 

ETJ^enia Patterson is a new senior clerk in the Publications Office* 

Carol Small has returned to join the Drafting Department o 

Terminated Non-Staff 

Lillian tiurso Nancy Heselton 

Eleanore Galant Gharlebert Ingram 

Sheila Heffeman Edmund Landers 



